Asymptotic properties of self-energy coefficients 
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We investigate the asymptotic properties of higher-order binding corrections to the one-loop self- 
energy of excited states in atomic hydrogen. We evaluate the historically problematic A§o coefficient 
for all P states with principal quantum numbers n < 7 and D states with n < 8 and find that a 
satisfactory representation of the n-dependence of the coefficients requires a three-parameter fit. 
For the high-energy contribution to ^60, we find exact formulas. The results obtained are relevant 
for the interpretation of high-precision laser spectrocopic measurements. 
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Bound-state quantum electrodynamics (QED) occu- 
pies a unique position in theoretical physics in that it 
combines all conceptual intricacies of modern quantum 
field theories, augmented by the peculiarities of bound 
states, with the experimental possibilities of ultra-high 
resolution laser spectroscopy. Calculations in this area 
have a long history, and the current status of theoret- 
ical predictions is the result of continuous effort. The 
purpose of this Letter is twofold: first, to present im- 
proved evaluations of higher-order binding corrections to 
the bound-state self-energy for a large number of atomic 
states, including highly excited states with a principal 
quantum number as high as n — 8, and second, to ana- 
lyze the asymptotic dependence of the analytic results on 
the bound-state quantum numbers. Highly-excited states 
(e.g, with n — 4 to 12) are of particular importance for 
high-precision spectroscopy experiments in hydrogen (for 
a summary, see for instance |l|, p. 371]). 




FIG. 1: Plot of the radial probability density r 2 |-0(r, 6», 0)| 2 
of the nonrelativistic 8D wave function (angular momentum 
projection m = 0) in the plane of constant azimuth </> — 
0. The calculation of the relativistic Bethe logarithm 
starts from this nonrelativistic wave function, with relativistic 
corrections being taken into account via generalized Foldy- 
Wouthuysen transformations 0, [?| . 

In the analytic calculations, we focus on a specific 
higher-order binding correction, known as the co- 



efficient or "relativistic Bethe logarithm." We write the 
(real part of the) one-loop self-energy shift of an electron 
in the field of a nucleus of charge number Z as 



AE SE = 



a (Za) 4 



F(nlj, Za) m c 2 



(1) 



where F(nlj, Za) is a dimensionless quantity. In this 
Letter, we use natural units with h = c = m = 1 and 
e 2 = Ana (m is the electron mass). The notation nlj is 
inspired by the usual spectroscopic nomenclature: n is 
the level number, j is the total angular momentum and 
I is the orbital angular momentum. 

The semi-analytic expansion of F(nlj, Za) about 
Za = for a general atomic state with quantum numbers 
n, I > 1 and j gives rise to the expression, 



F(nlj,Za) = A i0 {nlj) + [Zaf [A 61 {nlj) hx(Za)' 
+G SE {nl J ,Za)} (J>0), 



(2) 



where GgE {nlj , Za) — > constant as Za — > 0. The limit 
as Za — * of Gssfolj, Za) is referred to as the Aqq 
coefficient, i.e. 

.4, 



60 (nlj) = lim G S E{nlj,Za) . 

Za. — >0 



(3) 



It is this coefficient which has proven to be by far the 
most difficult to evaluate 0, E3 [1 IE El Furthermore, 
the complexity of the calculation increases sharply with 
increasing principal quantum number n, both due to the 
more involved structure of the bound-state wave function 
(see also Fig. QJ, and due to the necessity of subtract- 
ing bound-state poles that lie infinitesimally close to the 
photon integration contour. The atomic states with the 
highest n for which analytic results are available today 
are the 4P states 0. In this Letter, we present analytic 
data for the coefficient of P states with n < 7 and 
all D states with n < 8. For a given n, the calculation is 
more involved for nP than for nD, because there is one 
more term in the nonrelativistic radial nP wave function 
than in the corresponding nD wave function (when they 
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are expressed as a function of the electron-nucleus dis- 
tance). Essentially, the number of terms in the radial 
wave function determines the complexity of the calcula- 
tion. 

One of the most demanding specific calculations in the 
evaluation of A 60 is necessitated by a Bethe-logarithm 
type contribution given by the relativistic wave-function 
correction Fs^ this contribution is defined in Eqs. (43) 
and (53) of [8j. For 7P and 8D states, we use up to 
200,000 terms in intermediate steps in the evaluation of 
this correction. Because involves relativistic correc- 
tions to the coefficient A^q , which in turn is given mainly 
by the Bethe logarithm, it is natural to refer to the entire 
coefficient as a "relativistic Bethe logarithm." 

The "normal Bethe logarithm" \nko(nl) forms part of 
the coefficient A40 for which a well-known general for- 
mula (see, e.g., Ref. Q, p. 468]) reads 



A i0 (nlj) = - 



1 



2k, (21 + i; 



- - In k (nl) 



(4) 



where k = 2 (I — j) (j + 1/2). Formulas for A§\ valid for 
P and D states read as follows (see, e.g., Ref. 0, p. 468]) 



1 ( 2Q 
A 61 (nP 1/2 ) = - (33-^ 



Aei{nP 3 / 2 ) = |(9-^) 



(5) 



(6) 



A e i(nlj) 



32 3 - 



/(/ + !) 



3 n (2l + m) 

m— — l 



(l>2). (7) 



Note that Aei(nlj) — > constant for n — > 00 at constant I 
and j. It is the purpose of this Letter to present new re- 
sults for the Aqq coefficients. Details of our calculations 
will be presented in a forthcoming article llOl . It has 
been observed previously by Karshcnboim |ll| that the 
n-dependence of the Aqo(uP) coefficients can be fitted 
to a satisfactory accuracy by an (n 2 — l)/n 2 -type model, 
and a two-parameter fit has been employed for the n- 
dependence of the S-state coefficients A 60 (nS 1 / 2 ) [l2| . 
Our data for P states in Tab. |U are roughly consistent 
with this (n 2 — l)/n 2 model. 

For the atomic states under investigation, the self- 
energy contribution due to hard virtualphotons (high- 
energy part) obtained by the e method |6J, |g, 111 llfj is 



F H (nlj,Za) = 
+ (Za) 2 



1 



2k (21 + 1) 

JC Aei ln(2e) + 0(e) 

e 



(8) 



The ellipsis denotes higher-order terms, which are irrel- 
evant for the current investigation. In Eq. JSJ), /C and 



C, as well as Aq%, are state-dependent coefficients. For 
concrete evaluations of the high-energy part concernin 
specific atomic states, see H Eqs. (18) and (19)] and 
Eqs. (55)-(58)]. The low-energy part assumes the form 



F L (nlj,Za) 



(Za) 



-- In k (nl) 



£ H h Aei In , 

e \ (Za) 1 



(9) 



0(e) 



where we omit terms that are irrelevant at relative or- 
der (Za) 2 in the evaluation of F(nlj, Za). A detailed 
explanation of the e method will be given in [Io|. The 
dependence on C cancels when the high- and low-energy 
parts are added. Specifically, we have 

A 6 o =fC~ A 61 ln2 + £. (10) 

Upon inspection of JSJl and (J5J, we identify 

A 60 .„ =IC-A 61 In 2 (11) 

as the high-energy contribution to ^60 j and ^60,l = £ 
as the low-energy contribution (see Figs. and EJ. 
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FIG. 2: The plots show the dependence on the principal 
quantum number n of the high- and low-energy part of the Aeo 
self-energy coefficient, as well as their sum (^60)- The curves 
for the high-energy contribution Aao,H represent the exact 
results ©-10! I fl 11 ■ and ill 2t . with n being generalized to 
a continuous variable (only integer n values have physical 
significance). The smooth curves for the low-energy parts 
Aeo,L = C result from a three-parameter fit of the data in 
Tab.Qto the function in (1131 : the fit parameters are given in 
Tab.|H] The curves in the lower row represent the total result 
for ^60 = A 60jH + Ae,o,L- 
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TABLE I: A 60 coefficients for P 1/2 


, P3/2, D3/2 and D5/2 states 


n = 2, ... ,7 for P states, and n = 3, . 


. , 8 for D). The quantity 


C L 


5 implicitly defined in Eq. @ and represents the low-energy 


contribution to ^60 ■ 




n 


A 6 o(nP 1/2 ) 


£(nP 1/2 ) 


Aeo(nP 3 / 2 ) 


C{nP 3/2 ) 


o 

L 


-0.998 904 402(f) 


-0.795 649 812(1) 


— U.0O0 3(3 4b4(l) 


n ko 1 r 1 {-' *7on/i *\ 
— U.084 OlD /O0(lj 


Q 
O 


-1.148 189 956(1)° 


-0.944 288 447(1) 


—0 ^97 ^69 388H \ 


—0 693 ^Sfifi 497H 1 


/I 

*T 


-1.195 688 142(1) 


-0.997 810 211(1) 


(S^fl QA^ 7Qf\(~\ \ 
u.uou y^o / yu^i ) 


7*30 ^7Q 1 17f~\ \ 


c 



-1.216 224 512(1) 


-1.023 991 781(1) 


(SzL7 01 3 ^OOfl \ 


7zL7 fil ^ fi^fl \ 

U. / *± 1 u±o UOO^l ) 


u 


-1.226 702 391(1) 


-1.039 079 399(1) 


O 1 ^/l SQ3f1 A 

— u.ooo 104 oyo^ij 


O 7^ RQ7 AQQflS 


7 


-1.232 715 957(1) 


-1.048 800 134(1) 


fSfi9 097 Rfifin \ 


7fi9 f!99 Q^fifl \ 
u.i uz uzz you ) 


n 


^6o(nD 3/2 ) 


C{nD 3/2 ) 


^60(nD 5/ 2) 


£(nD 6/2 ) 


3 


0.005 551 575(1) 


0.021 250 354(1) 


0.027 609 989(f) 


0.019 188 397(1) 


4 


0.005 585 985(1) 


0.022 882 528(1) 


0.031 411 862(1) 


0.020 710 720(1) 


5 


0.006 152 175(1) 


0.023 759 683(1) 


0.033 077 571(1) 


0.021 511 798(1) 


6 


0.006 749 745(1) 


0.024 294 690(1) 


0.033 908 493(1) 


0.021 975 925(1) 


7 


0.007 277 403(1) 


0.024 645 479(1) 


0.034 355 926(1) 


0.022 264 036(1) 


8 


0.007 723 850(1) 


0.024 886 986(1) 


0.034 607 492(1) 


0.022 452 259(1) 



a We take the opportunity to correct a computational error for this 
result as previously reported in Ref. Eq. (96)], where a value of 
-1.14768(1) was given. 



TABLE II: Coefficients L\ , L 2 and L3 that result from a least- 
squares fit of the n-dependence of our data for £ in Tab. □( sec 
also Figs. |21 and |3J. The value of L\ from this global fit should 
approximate the limit lim n ^oo C(nlj) in Eq. If 3H . although it 
is not necessarily the best estimate. 



state 


Li 


L 2 


U 


Pl/2 


-1.082 


0.0966 


0.950 


P3/2 


-0.775 


-0.0232 


0.811 


D3/2 


0.0264 


-0.00952 


-0.0175 


D5/2 


0.0235 


-0.00568 


-0.0220 



We obtain the following general formulas for IC: 
637 1 767 



/C(nP 1/2 ) 

£(nP 3 /a) 
^(nD 3/2 ) 

AC(nD 5 / 2 ) 



1800 

2683 
7200 4 

157 



1 

4n 

1 

16n 



5400 n 2 
2147 



30240 
379 



18900 



3 

~ 80n 
1 

6to ~ 



5400 n 2 ' 
3007 



37800 n 2 

1759 
18900 n 2 ' 



(12a) 
(12b) 
(12c) 
(12d) 



All of these formulas are consistent with a limit IC —> 
constant as n — > 00 for constant I and j. The n- 
dependence of the nonrelativistic C(nlj) contributions as 
listed in Tab. [fl can be approximated very well using a 
three-parameter fit inspired by the above structure found 
for the high-energy IC contributions. We find 



C{nl j )*L l {l j ) + 



+ 



(13) 



where L\, L 2 , and L3 assume values as listed in 
Tab. [H] for the series of states under investigation. 



The n-dependence of the low-energy contributions C is 
smoother than the corresponding curves for the high- 
energy part (see Figs. [21 and 0). The excellent agreement 
of the fits with the numerical values of j460,l, together 
with our exact results for the high-energy part as given 
by Eqs. I©-©, (fTT|) . and (JT2J, could suggest a constant 
as n 
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FIG. 3: The analog of Fig.EJfor D 3/2 and D 5/2 states. The 
minimum in Aeo,H {riD 3 / 2 ) near n = 5 is determined by the 
exact formulas JJJ and I12cll . 



For Rydberg states with the highest I possible for given 
n (i.e., I = I = n — 1), our results are consistent with 



lim Aeo(nlj) = for I — n — 1, j 



l±l/2, (14) 
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which is plausible to suggest as a conjecture. The conjec- 
ture is indicated by the trend in the numbers (j = 1 — 1/2) 

• 74 6 o(3D 3/2 ) = 0.005 551 573(1), 

• A 60 (4F 5/2 ) =0.002 326 988(1), 

• A 60 (5G 7/2 ) = 0.000 814 415(1), 
as well as in the results (j = I + 1/2) 

• A 60 (3D 5/2 ) = 0.027 609 989(1), 

• A 60 (4F 7/2 ) = 0.007 074 961(1), 

• A 60 (5G 9/2 ) = 0.002 412 929(1). 

The magnitude of A§$(nlj) appears to decrease faster 
than 1/n. In general, relativistic corrections acquire 
at least one more inverse power of n when n = I + 1, 
j = n — 1 ± 1/2, and n large, than S or P states of the 
same n. This can for example be seen in the relativis- 
tic correction of order (Za) 4 to the Schrodinger-Coulomb 
electron energy [Eq. (2-87) of [l^]. 



E n j = m - 



(Za) 2 m 
2n 2 



{Za) 4 m 



2j + l 



_3_ 

8^ 



0[{Zaf 



For j = n — 1 ± 1/2, this relativistic term acquires an 
additional inverse power of n. Our results suggest that 
analogous statements hold for radiative corrections given 
by relativistic Bethe logarithms. 

We have presented results of a calculation of higher- 
order binding corrections to the one-loop self-energy for 
highly excited hydrogenic atomic levels (see Tab.HJ). Cal- 
culational difficulties induced by the more complex ana- 
lytic structure of the wave functions have been a severe 
obstacle for evaluations of relativistic Bethe logarithms 
at high n, and no prior results are available for Aqo for 
any state with n > 4 (see Ref. 0). Intermediate expres- 
sions contained up to 200,000 terms; without a computer, 
this work would have been impractical. Our calculation 
is split into a high- and a low-energy part. We find that 
the dependence of the low-energy contribution to on 
the principal quantum number of the atomic state un- 
der investigation can in many cases be represented accu- 
rately using a three-parameter fit [see Eq. Ijl3(l and the 
data in Tab. ITT]. As suggested by the exact formulas for 
the high-energy part given in Eq. i|12fl and the curves 
in Figs. |2] and |21 a fit with less than three parameters 
cannot be assumed to lead to a satisfactory representa- 
tion of A eo . Our final results for A eo are given in Tab. [I] 
We establish that the magnitude of decreases rapidly 
for Rydberg states with the highest possible angular mo- 
mentum for each principal quantum number. Our cal- 
culations improve the knowledge of the self-energy of an 



electron bound to a nucleus [see Eqs. CEJ-© and Tab.p. 
They are motivated by the dramatic ally increasing preci- 
sion of laser spectrocopy [lj, |l3> llil UM i which is rapidly 
approaching the 1 Hz level of accuracy. For the deter- 
mination of fundamental constants from high-resolution 
spectroscopy, frequency measurements of at least two 
different transitions have to be performed. Highly ex- 
cited, slowly decaying D states are attractive because 
they can be excited out of S states via two-photon reso- 
nance 

The authors would like to acknowledge helpful discus- 
sions with K. Pachucki, S. Jentschura and J. Sims. U.D.J, 
acknowledges support from the Deutscher Akademischer 
Austauschdienst (DAAD). E.O.L. acknowledges sup- 
port from a Lavoisier fellowship of the French Min- 
istry of Foreign Affairs, support by NIST, and a grant 
of computer time at the CINES (Montpellier, France). 
G.S. acknowledges support from BMBF, DFG and from 
GSI. The Kastler Brossel laboratory is Unite Mixte de 
Recherche 8552 of the CNRS. 



[1] 

m 

[3] 

[4] 
[5] 
[6] 
[7] 

[8] 

[9] 
[10] 

[11] 
[12] 

[13] 
[14] 

[15] 

[16] 

[17] 



P. J. Mohr and B. N. Taylor, Rev. Mod. Phys. 72, 351 
(2000). 

G. W. Erickson and D. R. Yennie, Ann. Phys. (N. Y.) 
35, 271 (1965). 

G. W. Erickson and D. R. Yennie, Ann. Phys. (N. Y.) 
35, 447 (1965). 

G. W. Erickson, Phys. Rev. Lett. 27, 780 (1971). 

J. Sapirstein, Phys. Rev. Lett. 47, 1723 (1981). 

K. Pachucki, Ann. Phys. (N. Y.) 226, 1 (1993). 

U. D. Jentschura, P. J. Mohr, and G. Soff, Phys. Rev. 

Lett. 82, 53 (1999). 

U. Jentschura and K. Pachucki, Phys. Rev. A 54, 1853 
(1996). 

U. D. Jentschura, G. Son", and P. J. Mohr, Phys. Rev. A 
56, 1739 (1997). 

E.-O. Le Bigot, U. D. Jentschura, P. J. Mohr, P. Indeli- 
cato, and G. Soff, The Bound-Electron Perturbative Self 
Energy of non-S States, submitted (2003). 
S. G. Karshenboim, private communication (1997). 
S. G. Karshenboim, Z. Phys. D 39, 109 (1997). 

C. Itzykson and J. B. Zuber, Quantum Field Theory 
(McGraw-Hill, New York, NY, 1980). 

B. de Beauvoir, F. Nez, L. Julien, B. Cagnac, F. Biraben, 

D. Touahri, L. Hilico, O. Acef, A. Clairon, and J. J. 
Zondy, Phys. Rev. Lett. 78, 440 (1997). 

C. Schwob, L. Jozefowski, B. de Beauvoir, L. Hilico, F. 
Nez, L. Julien, F. Biraben, O. Acef, A. Clairon, Phys. 
Rev. Lett. 82, 4960 (1999). 

J. Reichert, M. Niering, R. Holzwarth, M. Weitz, T. 
Udem, and T. W. Hansen, Phys. Rev. Lett. 84, 3232 
(2000). 

M. Niering, R. Holzwarth, J. Reichert, P. Pokasov, 
T. Udem, M. Weitz, T. W. Hansen, P. Lemonde, G. 
Santarelli, M. Abgrall, P. Laurent, C. Salomon, and A. 
Clairon, Phys. Rev. Lett. 84, 5496 (2000). 



